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We consider spin and charge flow in normal metals. We employ the Keldysh formalism to find 
transport equations in the presence of spin-orbit interaction, interaction with magnetic impurities, 
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I. INTRODUCTION 

Spin flow in nanostructures has recently attracted con- 
siderable interest in the scientific community 1 ' 2 ' 3 . The 
vision of magnetoelectronics and spintronics is to inject, 
manipulate and detect spins in nanostructures which can 
give new functionality in electronic devices. The spin flow 
can be controlled by e.g. external electric or magnetic 
fields. 

A Hall voltage builds up perpendicularly to the current 
flow under an applied magnetic field in normal metals 
due to the Lorentz force. The Hall voltage increases with 
applied magnetic field. Magnetoelectronic circuits are of- 
ten realized by using ferromagnets that can spin polarize 
the current flow. In ferromagnets, there is an anomalous 
Hall voltage proportional to the magnetization, e.g. a 
transverse charge potential, even in the absence of an ap- 
plied magnetic field. The anomalous Hall effect is caused 
by the spin-orbit interaction, which correlates the mo- 
mentum of the electron with its spin. This causes an 
dependence of the electron flow with the relative angle 
between its direction and the non-zero magnetic order 
parameter in ferromagnets In ferromagnetic metals, 
spin-orbit interaction is also a source of crystalline mag- 
netic anisotropy energies since the spin-orbit interaction 
couples the magnetization with the crystal structure. 

The scattering mechanisms responsible for the anoma- 
lous Hall effect are skew scattering^ and side-jump 
scattering 8 as well as the anomalous velocity operator 
due to spin-orbit interaction and impurity scattering. A 
schematic picture of the skew scattering mechanism is 
shown in Fig. ^ After scattering off the impurity po- 
tential, there is a spin-dependent probability difference, 
represented by small angles, of the electron trajectories. 
This leads e.g. to a slightly larger chance that electrons 
with spin up moves upwards and electrons with spin down 
moves downwards after scattering. The side-jump mech- 
anism is also caused by the combined spin-orbit and im- 
purity scattering, see Fig. [21 After scattering off the im- 
purity, a small "side-jump" , develops between the trajec- 
tories of electrons with spin up and down far away from 
the scattering center. Additionally, the spin-orbit inter- 
action does not commute with the electron momentum 
operator. This leads to an anomalous velocity operator 




FIG. 1: Schematic picture of the skew scattering mechanism. 
An incident electron with spin up (down) scatters preferrably 
with a postive (negative) angle. 
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FIG. 2: Schematic picture of the side-jump mechanism. The 
trajectory of the outgoing electrons is shifted to the upper 
(lower) side the scattering center at large distances for spin 
up (down) states. 

that can contribute to the Hall effects, as well. 

Spin-orbit scattering is also important in normal met- 
als. It is well known that it causes a loss of spin co- 
herence. Hirsch predicted the existence of a novel spin 
Hall effect 5 analogues to the anomalous Hall effect in 
ferromagnets and developed a phcnomcnological theory 
for the effect. In the absence of spin-orbit scattering, 
electrons with spin-up and spin-down scatter equally on 
non-magnetic impurities. However, as seen above, for 
nonzero spin-orbit interaction, when a current passes 
through the sample an imbalance between left-moving 
and right-moving particles is established, and an accom- 
panying transverse spin accumulation potential builds up 
in the system. 

Zhang computed this spin potential in the diffu- 
sive transport regime, and found that it should be 
measurable^. He considered longitudinal transport in 
a thin normal metal film, and computed the resulting 
transverse spin Hall voltage. In this regime, the spin- 
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accumulation is weak and he found that the spin Hall 
voltage is governed by the anomalous velocity operator. 
By using the framework developed by Hirsch and Zhang, 
effects of contact resistances on the spin Hall effect have 
also been considerecUS. 

Spin-orbit scattering is also important in n- 
and p-doped semiconductors, where Rashba-type 
spin-orbit couplingii leads to interesting spin Hall 

effeCtO 12 ! 13 ! 14 ! 15 ! 16 ! 17 ! 18 ! 19 ! 20 ! 21 ! 22 ! 23 ! 24 ! 25 ! 26 . The ffiOSt 

interesting case is spin transport in a two-dimensional 
electron system (2DES). The study of this effect is 
controversial. The debato 16 i 17 i 19 i 20 i 21 i 22 i 24 i 27 is focused 
on whether or not the spin Hall conductivity has a 
universal value <t s h = e/(8irh) and in what regime this 
result is applicable. Our study is complementary to 
these studies of spin Hall effects in semiconductors. 
In our normal metal case, the extrinsio 2 ^ spin Hall 
effects arise due to the spin-orbit at impurities. In 
semiconductors, spin-orbit interaction can be important 
even in ballistic systems 2 ^, in the absence of impurities, 
for systems with broken spatial inversion symmetry. 

In this work, we derive the spin Hall effects in the pre- 
sense of spin diffusion from the Eilenberger equation 30 
in presence of spin-orbit coupling 3 ^, magnetic impurities 
and non-magnetic impurities. We use the Keldysh Green 
function technique in the quasiclassical approximation. 
Our calculations go beyond the assumptions in Ref. 0, 
which only included effects of the anomalous velocity op- 
erator, in that we rigorously derive diffusion equations 
that also include effects of skew scattering and side-jump 
scattering. Our results agree with the results by Zhang in 
the limit he considered, where skew scattering and side- 
jump scattering can be disregarded. We also consider 
spin and charge transport in normal metals in another 
transport regime, phcnomcnologically treated by Hirsch, 
when normal metals are biased by ferromagnets, in which 
spins accumulate in the normal metal even in the absence 
of the spin-orbit interaction. We demonstrate that in this 
regime, skew scattering and side-jump scattering cannot 
be disregarded, and compute the spin Hall and spin-orbit 
induced charge Hall effects. In our analysis, we also con- 
sider transport through ferromagnet-normal metal inter- 
faces beyond the assumptions in Ref. ^3 by using the 
boundary conditions obtained within magneto-electronic 
circuit theory, which is necessary in dealing with non- 
collinear spin and magnetization directions. Finally, our 
formalism also explicitly incorporates the effects of scat- 
tering off magnetic impurities, which can reduce the spin 
Hall effects. 

Our paper is organized in the following way. In the 
next section, Section [H] we outline our starting point 
microscopic Hamiltonian, and explain the diffusion equa- 
tions for spin and charge flow that we have obtained 
rigorously by using the non-equilibrium Keldysh Green 
function approach in the quasi-classical approximation. 
We compute from the diffusion equations the spin Hall 
effects in some relevant experimental geometries in Sec- 
tion [ffi] Our derivation of the charge and spin diffusion 



equation is an important part of our work, and its details 
are given m Section |TV| Finally, we give our conclusions 
m Section El 



II. MODEL AND TRANSPORT EQUATIONS 

We consider quasi-particles in a normal metal that in- 
teract with non-magnetic and magnetic impurities and 
include the spin-orbit interaction. The Hamiltonian of 
the system is 



H = -T^V 2 + Mmp + ^so + V m , 

2m 



(1) 



where we below will introduce the various terms. Impu- 
rity scattering is represented by the short-range potential 



(2) 



where rj is the coordinate of the i -th impurity cen- 
ter and 7j is the strength of the scattering potential. 
It is assumed that the scatterers are disordered so that 
< V imp (r) >= and < V imp (r)V imp (r') >= 5{r-r')nj 2 , 
where n is the impurity density and 7 2 is the average fluc- 
tuation of the scattering strengths. 

The spin-orbit interaction is described by the Hamil- 
tonian, 



p 



h.c. 



(3) 



where a is the dimensionless spin-orbit coupling constant, 
kp is the Fermi wave number and h.c. denotes the her- 
mitian conjugate. 

Magnetic impurites are introduced by 



V sm = V sm (r)a-S{r), 



(4) 



where V snl (r) is the strength of the coupling of the itin- 
erant electron spin to the spin of the magnetic impurity 
S(r). 

We are interested in the transport properties of diffu- 
sive system where the system size is much larger than the 
mean free path. A rigorous method to obtain the correct 
diffusion equation is to start from a microscopic descrip- 
tion using the non-equilibrium Keldysh formalism. We 
consequently employ the Keldysh approach with two ap- 
proximations. First, we consider the quasiclassical ap- 
proximation, which is valid on length scales much larger 
than the Fermi wavelength, L> Xp. Second, we use the 
diffusion approximation which is valid when the system 
size is much larger than the mean free path, L 3> I = vft. 

The full derivation of the diffusion equation is an im- 
portant part of the present paper, but it is technically 
complicated and in order to make the paper more easily 
accessible we delay its derivation to the interested read- 
ers in Section llVI First, we show and explain the spin 
and charge diffusion equations as well as the expression 
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for the corresponding currents that we obtain. We intro- 
duce charge and spin- distributions jj, c and /x s , so that 
the charge density and spin-density are 

n(r) = iV /x c (r), 
s(r) = N ot i s (r), 

where N is the density of states. After considerably 
algebra outlined in Section Hvl we find that the resulting 
diffusion equation for the charge distribution functions is 
simply 



current: 



/'< 



= 0. 



(5) 



Similarly, the spin-distribution for small spin-orbit inter- 
actions, a -C 1, is governed by 



DV 2 ^ 



1 



1 



Ms 



(6) 



where D 



\v f t is the diffusion coefficient in terms of 
the Fermi velocity vf and the elastic scattering time r, 
r so is the spin-flip relaxation time due to the spin-orbit 
interaction, 



1 



8a 2 
9r 



(7) 



and T sm is the spin-flip relaxation time due to magnetic 
impurities, 



iim sm N S(S + l)v, 



2 

sm 



Here we have expressed the strength of the magnetic im- 
purity potential from Eq. Q in the momentum repre- 
sentation v sm , n srn is the concentration of the magnetic 
impurities, and S is the spin of the impurity. Thus, both 
the diffusion equations for charge and spin have the fa- 
miliar forms used exensively in the literature for spin and 
charge transport. 

Skew scattering, side-jump scattering and effects of the 
anomalous velocity operator are all contained in the ex- 
pressions for the current. We find that the total 2x2 
current in spin-space can be expressed as 



3 = 3o+ 3a 
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(9) 



where j a is the ordinary current without spin-orbit inter- 
action, ji av is the current due to the anomalous velocity 
operator, j ss is the current due to skew scattering and 
j s j is the current due to side-jump scattering. The full 
derivation of these currents is given in Section IIVI The 
charge and spin currents can be obtained by the trace 
of Eq. @ with the unit matrix and the Pauli matrices, 
respectively. In the limit of weak spin-orbit interaction 
(to lowest order in a), which is relevant for most normal 
metals, we compute the following contributions to the 



ej = -f- P V Mc + VQu s <x)] , 

ejav = o- — - V x fi s + er x V/x c , 
a 

e Jsj = cr-(erV)/x s , 
ejss = -aa&(Vfi s ) , 



(10) 

(11) 
(12) 

(13) 



where a ~ e 2 NoD is the conductivity. The anomalous 
current, j' av , contributes both to the spin and charge cur- 
rent. Contributions from skew scattering and side-jump 
scattering, j ss and j S j, only affect the spin-current. The 
current contribution j s j depends on the divergence of the 
spin accumulation in the direction of the current for any 
spin, which is the side-jump mechanism. The current 
contribution j ss arises from skew scattering. If we com- 
pare our results with Zhang^, we see that in addition 
to the contributions from the anomalous velocity oper- 
ator we include terms representing skew scattering and 
side-jump scattering. We demonstrate below that these 
additional contributes could correctly be disregarded in 
the geometry in Ref. 0, but that they the dominate spin 
Hall effects in other systems. We also derive expressions 
for the total current in case of arbitary a (see Section 



III. EXPERIMENTAL IMPLICATION. 

Let us now employ our theory to calculate the mag- 
nitude of the spin Hall effect in experimental relevant 
geometries. For simplify, we will consider the cases with 
small spin-orbit interaction (a <C 1), where the small a 
expressions for the current, (|TU|l . (|TT1) . l(T2*|l and (|T3|) . are 
valid. 



A. Thin metallic film. 

We consider first a pure normal metal as considered in 
Ref. |9j and shown in Fig. A thin film normal metal of 
length L and width d is attached with perfect contacts 
with zero resistance to a left reservoir with local chemical 
potential /il and a right reservoir with local chemical 
potential //r. In a pure normal metal system, there is 
no spin-accumulation in the limit a — > 0. That means 
that the spin-accumulation is small, being induced by the 
spin-orbit interaction. From (|1 1|) . (|12fl and (|13fl we thus 
see that contributions from the skew scattering and side- 
jump scattering to the current are of a higher order in 
the spin-orbit scattering than the anomalous current and 
can be disregarded. In the case of pure normal metals, 
in the weak spin-orbit interaction limit, the current can 
thus simply be expressed in terms of 



3o 



(14) 
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where the anomalous current simplifies to 



ejav « o- — -a- x V^c 



(15) 



The effect of the spin-orbit interaction is consequently to 
induce a transverse spin Hall potential. The magnitude 
of the spin Hall effect depends on the system size and 
geometry. If the system size is smaller than the spin- 
diffusion length, a spin accumulation cannot build up 
within the system, and consequently the spin Hall effect 
vanishes. Therefore, we consider the situation when the 
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FIG. 3: The thin metallic film with a contact to reservoir. 

system size is much larger than the spin-diffusion length, 
L 3> l s f = \/Dt s , where the total spin-flip relaxation 
time has contributions both due to spin-orbit scattering 
and magnetic impurity scattering, l/r s = l/r so + l/r sm . 
The solution of the diffusion equation is similar to the 
treatment in Ref. |9j. At distances larger than the spin- 
diffusion length from the reservoirs, the spin potential 
only depend on the transverse y-coordinate. In this 
regime, the general solution of the diffusion equation (JHJ 
has the form 



Ms 



c 3 e 



y/kt . 



-y/kt 



(16) 



where C3 and C4 are constants to be determined. We can 
determine the constants C3 and C4 from the boundary 
condition that there is no particle or spin flow across the 
transverse boundaries, e.g. the 2x2 current in spin-space 
must satisfy: 



j y (x,y = ±d/2) = 0. 



(17) 



Using the simplified equations (|14fl and l|15|) determining 
the current, we find C3 and C4 and the spatially dependent 
transverse spin potential: 



Formula (|18fl expresses the transverse spin Hall effect 
for thin metallic films with small spin-orbit interaction. 
There are accumulation of spins directed perpendicular 
to the film. As we see, the spin-Hall effect vanishes when 
scattering by magnetic impurities dominates the spin- 
orbit scattering, r sm <C t so , as expected. When magnetic 
impurity scattering is weak, and in the limit of a wider 
film than the spin-diffusion length d ^> l s {, we make the 
observation that the magnitude of spin Hall effect is " uni- 
versal" , ^ = 1$ (d/2) 



•M$(-d/2): 



u {A) 

/ sH 



3 fiR - Ml 
2 k F L 



(19) 



e.g. the spin Hall voltage does not depend on the spin- 
orbit interaction constant a. By "universal", we mean 
that the spin Hall potential does not depend on the 
strength of the spin-orbit scattering potential. This im- 
plies that, as long as the scattering off magnetic impuri- 
ties is weak, light metals (e.g Cu, Al) will give a similar 
spin Hall voltage as heavy metals (e.g Pt). Note that the 
spin Hall potential depend on the Fermi wave vector of 
the metal, kp, and thus weakly depends on the type of 
normal metal through this dependence. The reason for 
the "universality" is that although the spin-Hall potential 
is proportional to the spin-orbit scattering it only builds 
up within the spin-diffusion length which is inversely pro- 
portional to the spin-orbit scattering strength. 

This section also illustrates the differences between 
the extrinsic spin-orbit scattering off impurities in nor- 
mal metals, as treated here, and spin-orbit scattering in- 
duced by the Rashba term and impurities in the two- 
dimensional electron gas formed in semiconductor het- 
erostructures. In the two-dimensional electron gas, there 
is currently a controversy whether the spin Hall con- 
ductivity can reach a universal value in dirty systems, 
CsH = e/(8irh), independent of the spin-orbit scattering 
strength. In normal metals, we see from (|15|) that the 
spin Hall conductivity is not universal, but the spin Hall 
voltage (|19fl can be " universal" when scattering off mag- 
netic impurities is weak. 



Thin metallic film in tunneling contacts to 
ferromagnet and normal metal. 



m$G/) = 



3 M-R - HL 

8 k F L 



sinh(y/Z sf ) 



cosh(d/2/ s f) 



(18) 



This expression exactly agrees with the result of Zhang 9 
Ms,* = M T - M X = l s fE x C h /Csmh(y/l s f)/cosh.(d/2l si ) 
when we insert the Drude conductivity C = 
e 2 r(kp) 3 /6ir 2 m, the anomalous Hall conductivity Ch = 
e 2 a (kp) 3 /6n 2 (the dimensionless a = a hkp) and the 
electric field in the x-direction E x — (/j>r — /j,l)/L. In ad- 
dition, we use the identification between the spin-orbit 
spin-flip relaxation time and elastic scattering time, Eq. 
0- 



We will now consider another transport regime, in 
which there is a finite spin-accumulation present in the 
normal metal even in the absence of spin-orbit scattering. 
The spin-accumulation can be achieved by sandwiching 
the normal metal between two ferromagnets when the 
system is driven out of equilibrium, see Fig. There 
are two extreme regimes depending on the ratio between 
the system length and the spin-diffusion length. The first 
regime, when the length is large, (L ^> l s {), does not pro- 
vide novel physics since it is similar to the previous ex- 
ample. In this case any spins injected from the ferromag- 
nets will be lost, and the system resembles a pure normal 
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FIG. 4: The thin metallic film with tunneling contacts Ti and 
T2 to ferromagnet and normal metal. 



metal where a spin Hall voltage can build up transverse 
to the current direction. There is no charge Hall effect 
produced by the spin-orbit interaction in the middle of 
the system in this case. Therefore, we consider the sec- 
ond, more interesting regime, when the normal metal is 
short, (L -C l a f). In this case both spin and charge Hall 
potentials can build up transverse to the current direc- 
tion governed by the spin-orbit scattering. 

In order to further simplify the computation of the 
diffusion equations |J5}, JfjJ we assume that the normal 
metal is narrow, e.g. d <C l s { as well. The diffusion 
equations then simplify to V 2 /i c = and V 2 /x s = 0. 

In addition, we need the boundary condition of the spin 
and charge flow through the tunnel contacts from the fer- 
romagnetic reservoirs into the normal metal wire. These 
boundary conditions follow from magnetoelectronic cir- 
cuit theory, so that the interface transport can be de- 
scribed by spin-dependent conductances G^ and G^ for 
spin aligned and antialigned to the magnetization and 
a mixing conductance of reflection G^ for spins in the 
normal metal that are non-collinear to the magnetiza- 
tion direction^ For most systems, ReG^ 3> ImG^, and 
this will be assumed in the following, which simplifies the 
expressions for charge and spin currents^. The charge 
current along the transport, x, direction through the tun- 
neling contact can be written as 

el c = (G T + [ M f - /i c (0, y)] - (G T - G^m ■ Ms (0, y) 

in terms of the local quasi-equilibrium chemical potential 
in the ferromagnet [if and the spin and charge chemical 
potentials in the normal metal close to the ferromagnetic 
interface fi s (x — 0,y) and fi c (x — 0,y). Here m is the 
magnetization of the ferromagnet (| m | = 1). For a tunnel 
contact 2ReG n = G T + G^-M The spin-current along the 
transport, x, direction is then 

el s = m[(GT - G^f - /i c (0,y)) + (G* + G^f] - 
2ReG^fi s (0,y). (21) 

We assume that the spin-orbit interaction is weak, and 
expand the charge and spin accumulations to first order 
in a: 



^c{x,y) 
Vs(x,y) 



fi^{x,y)+aSfi c {x,y), (22) 
Hi°\x,y)+a5ti s {x,y), (23) 



where the zeroth order terms are the charge and spin ac- 
cumulations in the absence of spin-orbit interaction and 



the corrections are caused by the spin-orbit interaction. 
First, we compute the zeroth order terms that correspond 
to a = 0. The boundary condition jy{x,y = ±d/2) = 
of no current out of the transverse edges dictates that the 
general solutions /j,^ and only vary in the transport 
direction: ^ = c + cif , ^ = c a + c x f . 

Now we can equate the current through the tunnel con- 
tact with the current in normal metal at x = 0. This is 
a boundary condition on the ferromagnet T± (see Fig. 4) . 
For the charge current 

- G oCl = (G T + G x )( M l - c ) - (G T - G^m ■ c , (24) 

where the conductance of the normal metal is G = 
awd/L, w is the width of the film. For spin-current 

- G oCl = -(G T + G l )c + m[(G T - G l )(fi L - c ). (25) 

At the normal metal-normal metal tunnel contact Ti (see 
Fig. 4) we use the same expressions for currents (|2"Ull . (|2"TJ) 
with G = G^ = G^ , and the boundary conditions are 



G(c + ci - hr) 
G(c + ci) = 



= -G D ci, 
-G ci. 



After straightforward calculations we find 



^>(x) 



1 - 
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G + G L 



G 



^tot 



G 



G 



(fJ-i 



(26) 
(27) 



(28) 



/4 0) (z)=co 



Go \ Gl 



- r )( l iL-m) T , (29) 



where the totale spin-dependent conductances of the sys- 



tem are given by the resistances in series 1/G 



TCI) _ 

tot — 

l/GfiU' + 1/G + 1/G . Co can be expressed similarly, 
but is not shown here since it does not govern the spin- 
orbit induced correction to the currents which will be 
considered next. 

Next we introduce spin-orbit interaction which will 
produce current correction 6j(x,y) so that the full cur- 
rent is expressed as j = j° + Sj. The boundary condition 
$Jv(y = ±d/2) = gives 



5(i 



c \y=±d/2 



3mD dx s ' z 



dy 

2d 

g^ 6 ^,x \y=±d/2 

d d 



3dx d ^y 



dy 



y=±d/2 



dx 
h d 
3mD dx 



8nP 



(30) 
(31) 
(32) 
(33) 



where we have introduced the total polarization of the 
conductance of the system p to t = (Gj ot — Gt ot )/(Gj ot + 
Gj[ ot ). The derivative of rr-component of the correction 
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to the spin potential is governed by side-jump current, y- 
component by skew scattering current and z-component 
by anomalous current. In the limit d <C L we can expand 
the corrections to the chemical potentials in the small 
parameter (d/L), so that the solutions can be represented 
as 



non-zero only for nonzero magnetization of the ferromag- 
net in the z-direction (direction which is perpendicular 
to the film). At the same time the spin-Hall effect of the 
spins along z, l|41|l . is independent of the magnetization 
direction. So, assuming m = {0,0, 1} the magnitude of 
Hall effect is 



8Hc 
8 n s 



5^ - 



(1) 



(i) 



(34) 



d Gptot Mr - 
"iG + Go k F L 



(35) 

and the spin Hall effect 



(42) 



Due to simplified diffusion equations we write 



5 p. 



W _ M 

c — L 
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00 
c\ X 



00 



c\ x + c 2 y. 



(36) 
(37) 



(i) 

where z = 0,1. We are interesting for coefficients c 2 and 

because acj- 1 ' <C Cj/L, ac^ <§; cf, j — 0,1. The 
boundary condition on the contact 7\ gives 



i 6flc u=o 



dx 



|x=0 



G T + G L 

G u 
G T - G L 

G,, 
4 d 



5fj, c (0,y)m 



G T - G L 
G u 

gt + G l 



3 dx 



Su (0) e 

u Ms,ie c :c) 



and similarly for the contact T 2 . If we consider these 



conditions for b\i c = Sfj,c, S^i s = 5^' we derive c. 



(0) 



„(0) 



,(o) 



0. Conditions lj3U)l-{55fl give 



„(i) 
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c (1) 



fim z Gptot Ml - MR 
3mL> G + G d 
2 G ptot Ml - MR 



2m 
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Gptot Mi - MR 



c (1) - - 



G ~\~ G 
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G tot G 'tot Mi - MR 



3mI?GoG t T ot + Gi ot 



Finally, we have for the spin-orbit induced charge Hall 
effect 



5n[ B) (x,y) (xvn T 



Gptot Ml - MR 2/ 
G + Go fc F £ Z ' 



and for the spin Hall effect 
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(38) 

(39) 
(40) 

■V, (41) 



where I is the electron mean free path. As we see from the 
equation l|38|) the spin-orbit induced charge Hall effect is 



(B) 
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a— 



rA si I 

"~ r tot >J tot 



MR - ML 



1 Go GL 



kpL 



(43) 



The spin-orbit induced charge Hall effect vanishes when 
the polarization goes to zero, as should be expected. In 
comparison, we give the magnitude of spin Hall effect in 
pure normal metal regime in the limit of d <C l s f. As 
directly follows from eq. (|18fl 



M$(d«/ S f) 



a 



d fiR - Ml 
I kpL 



(44) 



To evaluate the expressions ()42[) . (|43|l and compare them 
to (|44|) we assume the conductance at the tunnel barrier 
T 2 equals to the sum of spin-dependent conductances at 
the barrier T\ , so G = G^ + G^ and take a typical value 
of the polarization of the conductance of the system p = 
(G T - G X )/(G T + G L ) = 1/2. The limit G > G D is not 
interesting because in this case our system will be similar 
to a pure normal metal attached to reservoirs, which was 
considered in the previous section. Another limit G <C 
G is also less interesting due to the small induced voltage 
across the normal metal and consequently vanishing Hall 
effects. So, we consider the most interesting case G = G . 
In this case 



^(d^l si ) 



(B) 



3 

7' 
1 

7 ' 



(45) 
(46) 



We see from these expressions that both the charge and 
the spin Hall effects are comparable to the spin Hall ef- 
fect in the pure normal metal regime attached to leads 
with perfect normal metal contacts in the regime L 3> l s f- 
Although the magnitudes of charge and spin Hall effects 
have a similar structure in both regimes, they have a dif- 
ferent origin in principle. The effects in this section are 
due to skew scattering, side-jump scattering and anoma- 
loys velocity, while the effects in the previous section are 
due to the anomalous velocity operator only. We treated 
in this section the regime d<SL L. We expect that the spin 
Hall and the charge Hall effects increase their magnitudes 

( B) (B) 

Mjj and Ms# with increasing of d until saturation when 
d ~ l s { similary to the case when the spin Hall effect in 
pure normal metal attains the "universal value" (|19|) . 
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IV. MICROSCOPIC DERIVATION OF 
DIFFUSION EQUATION 

We will in this section derive the diffusion equa- 
tions rigorously from the microscopoic Hamiltonian with 
the Keldysh Green's function technique in the quasi- 
classical limit. The nonmagnetic impurity potential 
(J2J) in terms of incident k and scattered k' wave vec- 
tors is Vi mp (fe, k') = 7^ exp — i (fe — fc') rj. Conse- 
quently, the spin-orbit interaction in this representation 
is V so {k, k') = -i(a/k 2 F )fr ■ (k x k') V imp (k - k'). Thus, 
in normal metals with dilute impurities the electrons in- 
teract with the potential V(k, k') = V5 mp + V so : 

V(k,k') = MV imp {k,k'), (47) 
where we have introduced the 2x2 matrix in spin-space 

(48) 



M = 1 -i- s -a{k x k') 



In addition, the electrons interact with magnetic impu- 
rities to be discussed below. Our transport theory is 
based on the Keldysh formalism 30,33 . In this formalism 
the Green's function has the form 



G 



G R G K 
G A 



(49) 



where the retarded, advanced and Keldysh Green's func- 
tions are 

G R = -W(t 1 -t 2 ){Mxi),9+{x 2 )] + ), (50) 
G A = +i9(t 2 -t 1 )(t*(x 1 ),y + (x 2 )} + ), (51) 
G K = -i{[*{xx),*+{x 2 )\J). (52) 

Here VP is the fermion annihilation operator, ^ + is the 
fcrmion creation operator, both in the Heisenberg pic- 
ture, and the coordinate Xi denotes both spatial position 
and time, Xj = (xi,£j). Note that the fermion annihila- 
tion and creation operators are 2-component vectors in 
spin-space. The self-energy has the same triangular ma- 
trix structure as the Green's function, 



£ = 



t A 



(53) 



We denote 4x4 matrices in Keldysh space by the sym- 
bol "check" and 2x2 matrices in spin space by the 
symbel "hat" Q. Next we define the center-of-mass and 
relative variables x = 3(xi + Xa), r = Xi — X2, and 
Fourier transform with respect to the relative coordinate 
r in order to obtain the Green's function in the mixed 
representation 

G(p,x) = y"dre- ip - r G(x + r/2,x-r/2). (54) 

We will also use the ^-integrated (quasiclassical) Green's 
function g(n,r) = (i/ir) J d£G(p,r), where £ = p 2 /2m — 



fi, and n is a unit vector along the momentum at 
the Fermi surface (n = kp/\kp\). Let us first con- 
sider the current produced by the normal velocity op- 
erator of the electrons, j'n- This current is expressed as 
Jn = ( e / 'm)di(^ + pty , where p is the momentum opera- 
tor. Introducing the Green's functions, the current due 
to the normal velocity operator is 

Jn(zi) = ~( hm (Vi - V 2 )G K (x u x 2 )) . (55) 

Inserting the Fourier representation of the Green's func- 
tion in the quasiclassical approximation gives for the or- 
dinary current in the mixed representation 3 ^ 



3n = J de J ^v F n (g K (n,r)) , 



(56) 



where (...) denotes averaging over impurities and N Q is 
the density of states at the Fermi level. 

The total current also has contributions caused by the 
anomalous velocity operator. This contribution to the 
current can be expressec 



Jav 



eN Q 



'I- I ^^(^o(n,n')g(n f ,r)) K . (57) 



The anomalous current caused by the spin-orbit interac- 
tion is 

v m (r) = = j^-t x W imp (r), (58) 

In Fourier space, 

v so (k,k') = -^r<r x (k | VMmp I k') = NV imp (k,k'). 
nkp 

where we have introduced the 2x2 matrix in spin-space 



JV = i-f-fr (As - fc') 

nkF 



(60) 



The challenge is now to find the average, (v so (n, n')g). 
Note that this is different than the average Green's func- 
tion appearing in the contribution from the normal veloc- 
ity operator (|5fc>|) . In the anomalous current (|58|l . a prod- 
uct of the spin-orbit scattering potential and the Green's 
function, both of which depends on the impurity config- 
uration has to be evaluated, and in general one should 
expect that (v so g) ^ (v so ) (g). Using the Dyson equation 

(<?) = .9o + .9oS (g) and (v so g) = NRI^ 1 (vg^ we find the 

result 



(v so g) = NM^H (g) 



(61) 



The self-energy part S can be expressed in the Born ap- 
proximation as S = (Vgv) and has one contribution 
due to scattering off non-magnetic impurties and two 
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contributions due to the spin-orbit interaction, S(n) = 
Ei(n) + S$(n) + U 2 o>(n) + £ sm (n), where 







(62) 


£(2) 

^SO 


2r Wnxn')j(n>(nxn'))„, 


(63) 


^so 


= ~ ^7 (9(n')a(n x n') + h.c.) n , . 


(64) 



Scattering by magnetic impurities (0} do not cause ad- 
ditional terms in the expressions for the current density, 
but gives an additional contribution the the electron self- 
energy 

S s ™ = -TT—l^ti (9(n')) n ,a u (65) 

2 

where the spin-flip relaxation time due to magnetic im- 
purity scattering is 

— = ^-n sm N S(S + 1) / ^-\V sm {n - n')\ 2 , (66) 

n sm is the concentration of the magnetic impurities and 
S is the quantum spin number of the impurity. The 



self-energy due to scattering off magnetic impurities will 
leads to simple additional spin-flip relaxation terms in 
the spin diffusion equation. In normal metals, the re- 
tarded component of Green's function equals g = 1, 
and g A = —g R = — 1. By evaluating the Keldysh com- 
ponent Eq. I|61|) the anomomalous current is given as 

Jav = Jde J £^NM-i[(2 + a 2 )g K (n') + 

a 2 &(n x n')g K (n')&(n x n') — ia{g K (n')a-(n x n') + 
a{n x n')g K {n')}\ + h.c. (67) 

We consider transport in the diffuse transport regime. In 
the diffusive regime, characterized by Vft <C L [L is the 
system size) the Green's function is almost isotropic 30 , 
and we can then expand the Green's function in form of 
isotropic and nonisotropic parts, 

J deg K (r,n)= fi a {r) + j\(r) n. (68) 

After integrating of the formula l|67[) we have the result 



Jav = { (a x 3i - 3i x 6)h +J2 ei Yl e ijk{Tr(vjh,k)f2 + Tr(a k j hj )f 3 ] ^ (69) 

jk 



fl = 


—a - 


a 

'T + 


h = 


2a 




T H 




h = 


2a 






-2a J 



where Tr(...) denotes the sum of diagonal elements of a matrix and the spin-orbit interaction strength dependent 
functions fi(a) are 

(n' z ) 2 -a(n z ) 2 (n' x ) 2 dn dn' 
1 + a 2 \n x n'\ 2 At: At: 
- 2aK) 2 (^.) 2 K) 2 -(n fc ) 2 (n;.) 2 dn^ 
1 + a 2 \n x n'| 2 At: At: 

2a{n l ) 2 {n' k ) i + 2a 2 {n k ) 2 {n' t ) 2 {n' k ) 2 - {n' k ) 2 dn dn' 
1 + a 2 \n x n'\ 2 At: At: 

I 



In the limit of weak spin-orbit interaction (a <C 1) we 
thus find the current contribution due to the anomalous 
veloicty operator 

J ' av = S(^ X "-" Xjl ) (?3) 

We have now found the full expression for the current 
density in terms of the distribution functions and we 
will now proceed to compute the diffusion equation. The 
starting point in calculating the Green's function is the 
equation of motion for the impurity averaged quasiclas- 



sical Green's function, the Eilenberger equation. 

hv F V(g)+i[±,(g)]_=0, (74) 

where vp is Fermy velosity. We will in the following omit 
the impurity average symbol and only implicity write av- 
eraging over momentum k or k'. We calculate Keldysh 
components of the various self-energy commutators. For 
scattering off non-magnetic impurities we obtain the well- 
known result: 

i[t u g] K =*(g K -{g K )) (75) 
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For spin-orbit scattering, which gives rise to side-jump 
and skew-skattering, using (&(n x n')cr(n x n')) , = 1 
we hnd 



Next, we find a second equation by averaging over dn the 
(|74|l multiplated by n 



and 





k ha 2 ' 




T 



+ (<x(n x n'^n'H™' x n)) 



VFTVflo + (1 + ^f-)jl 
f [{J\ X <x) + (<7 X £)' 

Eq. (|86|l can be solved to give 



(86) 



Similarly, using (cr(n x n')) = 0, we find £ ; 
Eio )A = and 



( 77 ) jl - -^F^[ i+2^/3 V ^° + ^lV(/Xs<r) - K 2 (&V)fl s 



where 



(AT 2 - K!)ct(V/x s )] 



v(i) 



K iha 



(g K (n')o-(n x n') - <5-(n x n')g K {n')) 

(78) 



In the diffusive transport regime, we use the representa- 
tion 



fi a = /i c l + /x s <t . 



(79) 



Employing i|76|) , we find after averaging over momentum, 

K 

where e is the 



the expression for Clio = i J de 
energy spectrum of the system: 



^so , g 



J 'so 



iha ' 
"37 
2ha 2 
3r 



n(ji x <r) + n(<7 x j ± ) , (80) 
/i s d- + (n<r)(n/x s ) + jin . (81) 



Similarly, for magnetic impurity scattering, from (|65|l we 
find 



yR = 

sm 



ih 
' Or 

— I S/ii 



1, E 



sm 



itl 
— ' sin 



sm o 
ZT< 



/'ft 



(/i c l - JLt s <x) 



and 



A n h . 

"sm = At s cr + Jin ■ 

Tsm *Tsm 



(83) 



(84) 



After substitution of l(7Efr. |STJ) 5 ||SDJ|, (Hj! into Eilen- 
berger equation 17411 and averaging over dn we find 

1 - /8a 2 1 \ 

^v F Vji + — + U s cr = . (85) 

I 



K t = 
K 2 = 



4a 2 



(3 + 2a 2 ) 2 - 4a 2 

2a(3 + 2a 2 ) 
(3 + 2a 2 ) 2 - 4a 2 



(87) 

(88) 
(89) 



We now use the representation l|68[) to find from i|56|) the 
current caused by the normal velocity operator 



Jn(t-) = — / (fi + jxn) = — - — n . (90) 



After substitution (|87(l into equation above we find finally 
that there are three contributions to the current density 
from the normal velocity operator: 

jxt(r) = jo(r) + jsj (r) + j ss (r) 



1 , (82) where 



Jo = 

Jsj — 
3ss — 



eN n D , 1 



2 L l + 2 a 2 /3 V ^ c 
27 + 36a 2 + 20a 4 



3(3 + 2a 2 )(l + 2a 2 )(3-2a 2 ) 
a(3 + 2a 2 )eN Q D 
(3 + 2a 2 ) 2 - 4a 2 

aeN a D 
~ 2a 2 + 2a + 1 



V(/* s <x)] (91) 

(d-V)Ma (92) 
<x(V/x 8 ) (93) 



Also using expression for ji ()87(l we can rewrite the for- 
mula for anomalous current (|69|l like 



iav(r) = -^{(^LV x^ + crx V Mc ) /i(a) + 

e* E,- fc (a)/ a (a) - tf 2 (a)/ 3 (a)) V fc /x Sij + (K 1 (a)f 3 (a) - K 2 (a)fr(a)) V jt i s ,k}} (94) 

I 



Diffusion equations for both spin and charge distribution functions can be derived directly from the equation 1851) 
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after substitution (|87() : 



V 2 Mc = 



(2a 2 + 3) 2 D 2 /x, 
(2a 2 + 3) 2 -4a 2 Ms r so 



4a 2 L> 



(2a 2 + 3) 2 - 4a 2 



(95) 



V(V Ms )- 



(96) 

In the limit of weak spin-orbit interaction (a <C 1) we 
obtain the simplified diffusion equations JHJ) and JHJ and 
the simplified expression for the current density © , (|10|l , 

on, jnj, and jnj. 



V. CONCLUSIONS 

We have derived diffusion equations for spin and charge 
flow in normal metals and the associated expression for 
the spin and charge currents. The total current con- 
sists of four contributions: Ordinary current, anomalous 
current, side-jump current and skew scattering current. 
These macroscopic diffusion equations, allows computa- 
tion of charge Hall voltages and spin Hall voltages in 
pure normal metals or hybrid ferromagnet-normal metal 
systems. 

We have considered two experimental relevant geome- 
tries and calculated Hall and spin-Hall voltages in the 
case of weak spin-orbit interaction. In pure normal met- 
als with no ferromagnetic contacts, there is no charge 



Hall effect due to the spin-orbit interaction, and the spin 
Hall effect is caused by the anomalous current, in agree- 
ment with the observation in Ref. 91 In this geome- 
try, we have made the additional observation that the 
spin Hall voltage reaches and "universal" value indepen- 
dent of the strength of the spin-orbit interaction, when 
spin-flip scattering due to spin-orbit scattering dominates 
spin-flip scattering due to magnetic impurities. When 
a spin-polarized current is injected into a normal metal 
film, both a Hall voltage and a spin Hall voltage exist. 
The magnitude of the Hall voltage is governed by side- 
jump, skew scattering and anomalous currents when the 
system is shorter than the spin-diffusion length. For sys- 
tems longer than the spin-diffusion length, the Hall volt- 
age vanishes, and the spin Hall effect resumes the value 
dominated by the anomalous current. In the interme- 
diate regime, the competition between skew scattering, 
side-jump scattering and anomalous velocity operator de- 
termines the spin Hall and charge Hall voltages. 
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